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Abstract 

This paper deals with the problem of precoding in multibeam satellite systems. In contrast to general 
multiuser multiple-input-multiple-output (MIMO) cellular schemes, multibeam satellite architectures 
suffer from different challenges. First, satellite communications standards embed more than one user in 
each frame in order to increase the channel coding gain. This leads to the different so-called multigroup 
multicast model, whose optimization requires computationally complex operations. Second, when the 
data traffic is generated by several Earth stations (gateways), the precoding matrix must be distributively 
computed and attain additional payload restrictions. Third, since the feedback channel is adverse (large 
delay and quantization errors), the precoding must be able to deal with such uncertainties. In order to 
solve the aforementioned problems, we propose a two-stage precoding design in order to both limit 
the multibeam interference and to enhance the intra-beam minimum user signal power (i.e. the one 
that dictates the rate allocation per beam). A robust version of the proposed precoder based on a first 
perturbation model is presented. This mechanism behaves well when the channel state information 
is corrupted. Furthermore, we propose a per beam user grouping mechanism together with its robust 
version in order to increase the precoding gain. Finally, a method for dealing with the multiple gateway 
architecture is presented, which offers high throughputs with a low inter-gateway communication. The 
conceived designs are evaluated in a close-to-real beam pattern and the latest broadband communication 
standard for satellite communications. 
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Index Terms 

Multibeam satellite systems, Precoding, Robust design. Multigroup multicast. 

I. Introduction 

A. Motivation 

Satellite eommunications will play a eentral role towards fulfilling next generation 5G eom- 
munication requirements m. As a matter of faet, anytime-anywhere conneetivity cannot be 
conceived without the presence of the satellite segment |l2l. Indeed, the satellite industry is not 
only targeting areas without backbone connectivity (maritime, aeronautic, rural), but also high 
dense populated scenarios with an existing communication infrastructure, where the satellite will 
become an essential element to decongest the terrestrial wireless network. 

In order to deliver broadband interactive data traffic, satellite payloads are currently implement¬ 
ing a multibeam radiation pattern. The use of a multibeam architecture brings several advantages 
in front of a single global beam transmission [l3]|. First, since an array fed reflector is employed, 
the antenna gain to noise ratio can be increased leading to high gain of each beam return link 
achievable throughput. Second, different symbols can be simultaneously sent to geographically 
separated areas, allowing a spatially multiplexed communication. Last but not least, the available 
bandwidth can be reused in sufficiently separated beams, leading to an increase of the user 
bandwidth yet maintaining a low multiuser interference. 

Nevertheless, whenever the system designers target the terabit satellite system (i.e. a satellite 
system offering a terabit per second capacity), the aforementioned multibeam architecture shall 
be reconsidered. Precisely, full frequency reuse among beams becomes mandatory in order to 
support the terabit capacity as larger available user bandwidth is required. As a consequence, 
when considering the satellite forward link, interference mitigation techniques need to be im¬ 
plemented either at the user terminal (multiuser detection) or in the transmitter (precoding). 

Whenever precoding is employed, the users must feed back their channel station information 
(CSI) to the transmitter so that it can revert the interference effect at the transmit side. These 
techniques relay severely on the quality of the CSI and they dramatically decrease their per¬ 
formance in case CSI is either deprecated or corrupted. On the contrary, multiuser detection 
techniques does not depend on the feedback channel but; however, the user terminal complexity 
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increases so as its eost. In addition, preeoding system level studies are providing eneouraging 
results for implementing this teehnique in real multibeam satellite systems [jH. As a result, we 
will eonsider preeoding as interferenee reliever for the present study. 

B. Previous Works 

The first designs of preeoding teehniques for multibeam satellite systems ean be found in [l5]|. 
Mimieking the linear techniques for multiuser multiple-input-multiple-output (MIMO) schemes 
the authors propose a zero foreing (ZF) and minimum mean square error (MMSE) preeoding 
designs. In addition, several ehallenges for the implementation of preeoding in multibeam satellite 
systems were pointed out. We deseribe them in the following. 

As a general statement, the payload eomplexity shall be maintained low and; eonsequently, 
the ground segment must perform most of the eomputations and transmit the preeoded signals 
through the feeder link. This feeder link must be able to support the overall satellite traffie, 
leading to a very large bandwidth requirement. This requirement is even larger when preeoding 
is deployed at the gateway sinee the feed signals must be preeoded and simultaneously transmitted 
through the feeder link. On the other hand, if the payload is equipped with a beamforming mech¬ 
anism, the feeder link bandwidth requirements ean be alleviated sinee only the user signals shall 
be transmittecji]. However, as presented in [|71 on board proeessing limits the overall gains 
obtained by the on ground proeessing [Ull, [[H. Remarkably, in case the on board processing is 
optimized eonsidering an underlying preeoding seheme, eertain throughput gain ean be preserved 
El. Yet another alternative is to shift the feeder link to the Q/V bands although certain diversity 
sehemes must be deployed for dealing with the large path loss EH- 

In order to deal with the extremely inerease of feeder link bandwidth requirements resulting 
from the full frequeney beam reuse pattern, several Earth stations (gateways) ean be deployed. 
With this, the available speetrum for the feeder link can be reused among spatially separated 
gateways through very direetive antennas. In these systems, the traffie is generated at different 
gateways so that several feeder links simultaneously transmit precoded data from isolated areas. 
As a consequenee, the preeoding teehnique shall be reeonsidered sinee, in order to have enough 

'This statement assumes that the number of feeds is larger than the number of users as it happens in most of the current 
deployments. This will be discussed in Section II. 
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spectrum to access all the feeds, each gateway has only access to a certain set of feeds. In 
[fT^ . lfT3ll a precoding scheme based on the regularized ZF scheme is presented. As a matter 
fact, the overall throughput is reduced when multiple gateways are considered even though 
computationally complex on ground schemes are deployed lfT4l . 

Last, but not least, satellite communications operate in a multicast fashion since data from 
different users is embedded in the same frame. Precisely, in order to increase the coding gain, 
each beam simultaneously serves more than one user by means of transmitting a single coded 
frame. This scheme entails a modification of the overall precoding scheme since more than one 
spatial signature per beam must be considered and; moreover, the achievable rate is dictated by 
the user with the lowest signal-to-noise-plus-interference ratio (SNIR). The first approach for 
solving this problem can be found in ifTSll were a sequential beamforming scheme is presented. 
However, whenever a large number of beams over the coverage area is targeted, one-shot scheme 
must be implemented. An example of this can be found in ifT^ and ifTTll where a design based on 
MMSE precoding and channel averaging is presented. Since the multicast multibeam transmission 
can be cast as multigroup multicast mathematical model, the proposed designs in ifT^ can be 
applied. Finally, very high throughputs can be obtained whenever the joint precoding and user 
scheduling is performed as in lfT9]l . 

C. Contributions 

In contrast to the aforementioned works, the present paper proposes a low complex ground pre¬ 
coding scheme to deal with the multibeam interference in multicast transmissions. The presented 
novel technique is based on two stage linear precoding similar to [|20]| . where the multiuser MIMO 
scenario is targeted. Our proposal offers higher spectral efficiencies than the regularized ZF [fT^ 
and the average MMSE scheme ifT^ . ifTTll yet maintaining a low computational complexity. 
It is important to remark that we prioritize low complex one-shot design in front of iterative 
gradient-based schemes lfT4ll . ifT^ . 

In addition, considering that the CSI will be corrupted at the gateway, a robust scheme is 
presented based on the first order perturbation theory of the eigenvectors and eigenvalues [|2T]| . 
This robust design is novel and it has not been applied before. The resulting precoding design 
remains low complex so that it can be implemented even if a very large number of feeds are 
considered. 
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Since the achievable rates deerease whenever the user channel veetors within one beam are not 
eollinear, we propose a user grouping teehnique. With this, over all possible users to be served 
for eaeh beam, we seleet the most adequate ones using a variation of the fc-means algorithm. 
This algorithm differs from the one presented in [fT^ as not only the ehannel magnitudes but 
the phase effeets are eonsidered. Additionally, a novel robust user grouping seheme is proposed 
in order to deal with the possible ehannel uneertainties. 

In ease the data traffie is generated by several gateways, a preeoding meehanism is presented 
for dealing with the main ehallenges; namely, CSI sharing and the distributed preeoding matrix 
eomputation. Both a redueed inter-gateway eommunieation for CSI sharing and a preeoding 
matrix division among gateways are presented. Even though the aehievable rates are deereased 
when the multiple gateway arehiteeture is eonsidered, the proposed seheme offers a good trade¬ 
off between eommunieation overhead, payload eomplexity and overall throughput. 

The rest of the paper is organized as follows. Seetion II introduees the problem so as the 
ehannel modelling. Based on this model and problem statement, Seetion III presents the two 
stage preeoding design for dealing with both the multibeam interference and intra-beam signal 
enhaneement. Relaying on a eertain preeoding design, a robust seheme is proposed based on the 
first order perturbation method in Seetion IV. Seetion V shows how to implement preeoding in a 
multiple gateway arehiteeture. Seetion VI presents numerieal simulations eonsidering the digital 
video broadband S2X (DVB-S2X) and a beam pattern of 245 beams. Seetion VII eoneludes. 

Notation: We adopt the notation of using lower ease boldfaee for veetors, v, and upper ease 
boldface for matriees, A. The transpose operator and the eonjugate transpose operator are denoted 
by the symbols (■)^, (•)^ respeetively. E[-] denotes expeetation. I denotes the identity matrix. C 
denotes the eomplex numbers. || ■ || denotes the Euelidean norm. | ■ | denotes the absolute value. 
A denotes the eomponentwise inequality, o denotes the Hadamard produet. 

II. Problem Eormulation 

A. Channel Model 

Eet us eonsider a multibeam broadband satellite with fixed reeeivers where a single gateway 
is provisioning signals to be transmitted through a feeder link. Over the feeder linlj^ a total 


^The feeder link connection is considered ideal. This is, noiseless and without channel impairments. 
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number of N feed signals are frequency multiplexed so that the payload has to detect and route 
them through an array fed reflector. This array fed reflector transforms the N feed signals into 
K transmitted signals (i.e. one signal per beam) to be radiated over the multibeam coverage 
area. 

As a matter of fact, the array fed reflector can have a single-feed-per-beam (SFPB) architecture 
whenever = iT or a multiple-feed-per-beam (MFPB) when N > K. This latter payload 
architecture presents lower beamforming scan losses and larger antenna gains than the SFPB 
[l22]| . For the sake of generality will consider the MFPB structure in the rest of the paper. 

Towards a spectrally efficient communication, all beams share the frequency band and; in a 
given time instant, the /c-th beam simultaneously serves a total amount of Qk users. In other 
words, in a give time instant the scheduler selects a set of Qk users at the /c-th beam (same for all 
beams) and it constructs a codeword with information to be transmitted to all Qk- Without loss 
of generally, we will assume that each beam serves the same number of users simultaneously 
and it is equal to Q (i.e. Qk = Q k = 1,, K). Figure [U shows the overall system. 



Fig. 1. The picture depicts the multicast multibeam satellite structure. The gateway delivers certain data to the coverage area 
by first the feeder link and; posteriorly, the satellite. While the feeder link multiplexes N signals, the satellite that is equipped 
with an array fed reflector, radiated a total of K signals (one signal per beam). Every radiated signal by the satellite shall be 
detected by a total number of Q users per beam. 


Under this context, the received signal can be modelled as 


y = Hx + n 


(1) 
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where y G is a vector containing the received signals at each user terminal. Vector 

n G contains the noise terms of each user terminal and we will assume that they are 

Gaussian distributed with zero mean, unit variance and uncorrelated with both the desired signal 
and the other users noise terms (i.e. E [nn^] = \kq)- 
The channel matrix can be described as follows 


H = F o (2) 

Matrix F G represents signal attenuation generated via both the atmospheric fading and 

the antenna feed radiation. This matrix can be decomposed as follows 


F = AG 


( 3 ) 


where A G is diagonal matrix whose diagonal entries are the atmospheric fading terms 

corresponding to the g-th user in the A:-th beam. Matrix G G takes into account the rest 

of gain and loss factors. Its {kq,n)-\h entry can be described as follows 


^ R^kqn 


- — 4 ^^^== ( 4 ) 

with dkqn the distance between the g-th user terminal in the k-th beam and the satellite. A is the 
carrier wavelength, is the Boltzmann constant, Bw is the carrier bandwidth, G\ the user 
terminal receive antenna gain, and Tr the receiver noise temperature. The term akqn refers to 
the gain from the n-th feed to the g-th user in the k-th beam. It is important to mention that the 
G matrix has been normalized to the receiver noise term. 

Furthermore, matrix $ G represents the phase variation effects between the n-th feed 

signal and the g-th user located in the A:-th beam. Each entry is defined as 


1 * 1 ,( 5 ) 

where 9g^n is uniformly distributed between 0 and In. 

The aforementioned phase effect determines the performance of the system severely. Due to 
that, an option is to use ultra stable oscillators so that the phase variation due to the different 
feed oscillators is reduced. With this, we can assume another phase distribution such that 

C ™ = ^<1 + ( 6 ) 


where (pq is a uniformly distributed in 0 and Itt and 7 ^ „ is modelled as a zero mean Gaussian 
distribution with variance y. 




For notation convenience, it is better to define G as the channel matrix for the k-th 

beam so that 

H= (7) 


Moreover, the channel vector of the g-th user in the k-th beam is defined as g so that 

Hfc= [hl„hl2,---,KQf ■ 


( 8 ) 


In order to minimize the multiuser interference resulting from the full frequency reuse, linear 
precoding is considered. Under this context, the transmitted symbol can be modelled as 


X = Ws (9) 

where s e is a vector that contains the transmitted symbols which we assume uncorrelated 
and unit norm (U = Ix) • Matrix W G is the linear precoding matrix to be designed. 


B. Precoding Design 


Let us formulate the precoding design of a multicast multibeam satellite system. The overall 
system performance can be optimized considering the maximum sum rate: 

K 

maximize > minimum „ 

W q=l,...Q 

k=l 

U- ( 10 ) 

subject to 


Tr (WW^) < Pt 

where ^ denotes the achievable rate of the g-th user at the k-th beam. 


rk,q = log2 1 + 


E 


K 
jjtk 




( 11 ) 




Pt denotes the available power at the satellite and corresponds to the k-th column of matrix 
W. It is important to remark that even though power sharing mechanisms among beams can be 
implemented [l2^ . [|2^ their deployment into next generation satellite payloads will require costly 
and complex radio-frequency designs. Under this context, the precoding design shall consider a 
per-feed available power constraint so that the constraint in (fTOl) becomes 

Pt 


(WW^) < 

\ / n.n — 


N 


n = 1, 




(12) 


where it has been assumed that the available power is equally share by all feed elements. 



In any case, problem (fTOl) is a diffieult non-convex problem whose convex relaxations even 
require eomputationally demanding operations GSll . For multibeam satellite systems the eompu- 
tational eomplexity of the preeoding design is an essential feature sinee these systems usually 
operate with hundreds of beams. Consequently, in contrast to other works, the target of this paper 
is to design a low eomputationally complex preeoding seheme able to achieve high throughput 
values. This is presented in the next seetion. 


III. Generalized Multicast Multibeam Precoding 

The preeoding design in multieast multibeam satellite systems has to main roles. First, the 
inter-beam interferenee shall be minimized and; seeond, the preeoding shall inerease the lowest 
SINK within each beam. Under this context, the precoding design ean be divided into two 
sub-matriees such as 

W = aW.Wfe, (13) 


where 


(14) 


and 


Wb 


Wfei 0 

0 Wb2 


0 

0 


e 


xK 


(15) 


_ 0 0 ■ ■ ■ 

with Wqj, G and G The matrix is used to mitigate the inter-beam 

interference first and then the matrix W;, is used to optimize the intra-frame data rate (i.e. 
the rate of the served users) so that and denotes the preeoding for optimizing the rate 
at A:-th beam. Finally, the parameter a is chosen to set the available power constraint (both for 
the per feed and total power eonstraint). 

In the following subsections, two different designs for and a single design for W;, are 
presented. 
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A. Inter-beam Interference Mitigation Precoding 

1) Multibeam Interference Mitigation (MBIM): Let us define as 

Hi 


H 


k — 


Hfc-1 

Hfc+i 


H 


K 


^ f^iK-l)QxN_ 


( 16 ) 


First we observe the interference impact of the precoding matrix Wq. This can be done by 
means of constructing the equivalent combined channel matrix after the precoding effect: 


(17) 


where the k-th beam effective channel is given by and the interference generated to the 

other users is determined by HfcWa^,. 

As described in [l20ll an efficient design of is given by the optimal matrix of the following 
modified MMSE objective function 



HiW,, 

HiW„, • 

HiW,^' 

HW,= 

H2W,, 








minimize E 

Wa 


K 


Y. 


k=l 


H-n 

' Pt " " 


(18) 


where the term ^ is obtained considering a total power constraint. 
The solution of this optimization problem is given by 


W 


MBIM 




where 


and 


= Vfc, 


D,, = ( Sfc + 


- 1/2 


(19) 


( 20 ) 


(21) 


Note that it has been considered the singular value decomposition of 


tH 
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2) Regularized Zero-Forcing: Let e £^kqxkq (jg^ote the regularized channel defined as 

( 22 ) 


H(«) = + ^1, 

Pt ’ 


where the same regularization factor as that of the multicast-aware regularized zero-forcing 
is considered. Let us define as 


H 


(R) 


Hi 


(R)- 


h(^) 

^k-1 

{R) 
k+1 


H 


H 


(R) 
K . 


^ (^K{Q—l)xKQ 


(23) 


The SVD decomposition of matrix can be described as 


H 


(R) 


= uf s™ 


^ k ’ ^ k 


1 H 


(24) 


where it is emphasized that there is always a null space of dimension Q spanned by This 


matrix is used for this scheme such as 




e C 


NxQ 


(25) 


B. Intra-beam Precoding 

After the first precoding scheme, W^, the /c-th beam observes an equivalent channel = 
HfcWa^ G . Based on Z^ the system designer shall construct w^j.. Let us mention that the 

equivalent channel for the g-th user located at the k-\h beam is denoted by G so that 

■ ■ ■ )■ (26) 


A suboptimal yet efficient approach is to maximize the average SNR considering the equivalent 
channel matrix. This is done with the following optimization problem 


Q 

maximize > I zF Wb. I ^ 

q=l 

(27) 

subject to 

llwbJI" < 1. 
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Since the objeetive funetion in (1271) ean be re-written as ||ZfcW6j^||^, the optimal solution of (ITT]) 
is given by the eigenveetor assoeiated to the largest eigenvalue of matrix Z^. Note that this 
preeoding design offers a low eomputational eomplexity sinee only the eigenveetor assoeiated 
to the largest eigenvalue needs to be eomputed. 


C. Power control 


Onee both preeoding schemes are eomputed, it is time to ealeulate the value of a in order 
to fulfil the maximum transmit power eonstraints. In ease the the maximum per feed available 
power eonstraint is eonsidered, 


Pt 

Oi , . , 

A^maxi ((WW^). .j 

whereas in ease total power eonstraints are eonsidered 

Oi — 

NTr ■ 


(28) 


(29) 


IV. Robust Multicast Multibeam Precoding 


As it will be shown in the simulation seetion, the MBIM sehemes offers larger aehievable 
rates than the R-ZF. Consequently, we will eonsider this design so as the average optimization 
intra-beam method. In any ease, note that preeoding performanee relies on an aeeurate CSI fed 
baek by the reeeiver. However, this CSI suffers from eertain degradation due to quantization and 
transmission delay. Due to that, it is eonvenient to reformulate the optimization problem in order 
to take into aeeount these possible variations GTll . This is done in the following subseetions for 
both the inter-beam and intra-beam preeoding. 

First, let us introduce the perturbation where the transmitter do not longer have access to H 
but to a degraded version sueh as 

H = H + A, (30) 


where A G is the perturbation matrix where it is assumed to be eonstrained so that 


|A||^ < 7 . 


(31) 


For notational eonvenienee, it is important to define the following sub-matriees 


A=[AlAl...,AlY , 


(32) 
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Ai 


(33) 


— 


Afc_i 

Afc+i 




g (^iK-l)QxN^ 


(34) 


where G is the perturbation associated to the fc-th beam, Sk,q G is the pertur¬ 

bation associated to the g-th user located at the A:-th beam. Under this context, the perturbation 
sub-matrices can be constrained as follows 


|AJ|2 


< 7fc, 


^k\\~ < 7fc 


for A: = 1 itT and where 


K 

Ik = X^7i- 

l^k 


(35) 

(36) 


(37) 


Additionally, a definition of a lower bound of the perturbation matrices A^ is convenient 


|Afc||^ > 7fc. 


(38) 


Remarkably, finding the adequate 7 bounds for all different perturbation matrices is not an 
easy task. Indeed, the computation of the different bounds shall be done on a empirical basis 
considering the different error sources and their final value on the perturbation matrix. This study 
is out of the scope of this paper and we will only provide a sensitivity analysis in the simulation 
section. 


A. Robust Inter-beam Precoding 

Considering the MBIM scheme presented in the previous section, whenever the robust worst- 
case optimization problem (fTSl) is targeted, the following new optimization problem shall be 
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considered 


minimize 


maximize 

K 

k=l 


■ K 

E 

.fc=i 


|H,W„J|2 + 


KQ 

Pt ' 


n 


subjeet to 

||Afc||^<^ k = l,...,K, 


(39) 


where = Hfc +A^. Considering that the optimal design on (fTSl) leads to the eomputation of an 
eigendeeomposition, the perturbation matrix will both impaet the eigenveetors and eigenvalues. 
Robust designs generally only eonsider the effeet on the eigenvalues [|^ Chapter 7]; however, 
the impaet on the eigenveetors eannot be eonsidered negligible IITTI . 

The following theorem provides an approximate solution of the optimization problem in (l39l) . 
Theorem 1 The optimal inter-beam precoding matrix which approximately minimizes (l39l) is 


where 


and 


where 


Maj. = Vfc 

.-. n-1/2 

Daj. = \^k + Cfclj 
(^k = yl + 27fcCrmax Hfc) • 


(40) 


(41) 


(42) 

(43) 


The rest of the matriees are defined in Appendix A and not ineluded here for the sake of brevity. 
Proof: See Appendix A. ■ 

Remarkably, whenever inereases, the resulting robust preeoding design is more different 

than the original design In any ease, the eomputational eomplexity of the robust design 

remains the same. 


B. Robust Intra-beam Precoding 

Similarly to the previous robust design the intra-beam preeoding shall eonsider tentative 
perturbations on their channel matriees. Considering the average optimization seheme, worst-ease 
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robust optimization for the /c-th beam ean be deseribed as the following optimization problem 

maximize minimize 1111 ^ 

Wbj. Afc 

subjeet to 



where 

= + AfcW„,. (45) 

The next theorem presents an approximate solution of the aforementioned problem. 

Theorem 2 An approximate solution of (l44l) is which is the first column vector of matrix 

Lfc {ukN o (Lf Lf Tfc + TfcLf Lfc) + l) (46) 

where 

(^k + 1) (47) 

and the rest of matrices are defined in Appendix B. 

Proof: See Appendix B. ■ 

Note that whenever Uk inereases, the solutions is more different than the original solution. 
Again, the eomputational complexity remains the same as the non-robust case. 

V. User Grouping 

One of the main limiting factors of multibeam multicast precoding is the spatial diversity of 
the users within each beam. Indeed, whenever the targeted users in each beam have orthogonal 
channel vectors, this is, 

hjmhfc.n = 0, (48) 

for the fc-th beam for m,n = 1... ,Q and m re, the intra-beam precoding is not able to deliver 
the intended symbols. Under this context it is beneficial that the system designer performs a user 
grouping before the precoding matrix is computed so that users with collinear channel vectors 
are simultaneously served 
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A. k-User Grouping 

In a given time instant, the seheduler determines a set of tentative users to be served. The 
number of these scheduled users is considered the same for each beam, fixed and equal to Q. 
For each beam, obtaining the most adequate groups of users is a cumbersome problem. Note 
that first, the system designer shall determine the adequate number of groups Gk per beam and, 
posteriorly, group them into those Gk groups. Clearly, the overall system throughput will depend 
on the user density over the coverage area: the larger number of users over the coverage, the 
larger throughputs can be obtained. 

In order to solve this problem, we will consider a random pre-processing. This pre-processing 
consists of randomly choose a user from the beam and, posteriorly, obtain the group of users. 
Note that with this first processing, we are severely levering the computational complexity of 
the technique. 

Under this context, let us consider that we have elected an arbitrary user m within the k-th 
beam. The user grouping scheme shall obtain the closest Q — I users in terms of Euclidean 
norm. Mathematically, 

minimize 

nGl,...,(3 

subject to (49) 

n ^ m. 

Since the considered Q in is not expected to be large, this optimization only requires a set 
of Euclidean norm comparisons between the scheduled users. As it happens with the precoding 
design, the user grouping scheme suffers from degradation whenever the user channel vectors 
are corrupted. A method for robustly group them is presented in the next subsection. 

B. Robust k-User Grouping 

Whenever the channel vectors are corrupted by a certain perturbation, a worst-case optimiza¬ 
tion shall be performed 
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minimize maximize \\h.k,m + 8k,m-'^k,n - 8k,n\\ 

subject to 

n ^ m, 

II ^fc,5 II ^ Tfc,g' 

Next theorem provides an approximate version of the aforementioned problem. 

Theorem 3 An optimization problem whose solutions upper bound the original worst-case 
robust grouping problem (l50l) is 


minimize ||hfc,m - hfc,nir + 7n,q 

subject to 

n ^ m. 

Proof: It is a simple derivation considering the Cauchy-Schwarz inequality and the fact that 
given a randomly chosen user m its perturbation does not influence the grouping optimization. 


With this optimization it is evident that whenever uncertainty is assumed in the channel vectors, 
this shall be considered in the user grouping design by means of an additional scalar penalty. 
Remarkably, in case we consider the same uncertainty to all users, the proposed approximate 
solution remains the same. 

VI. Multiple Gateway Architecture 

As a matter of fact, there might be the case where the feeder link cannot support the overall 
satellite data traffic. For instance, a payload equipped with N = 155 feeds with a user channel 
bandwidth of 500 MHz requires a feeder link bandwidth of 77.5 GHz which is an unaffordable 
requirement even if the feeder link carrier is located at the Q/V band. 

In order to solve this problem, the feeder link might benefit from certain spatial reuse so that 
several gateways can simultaneously send the data to be transmitted over the satellite coverage 
area. Under this context, G gateways can reuse the available bandwidth for the feeder link leading 
to large increase of the user bandwidth (see Figure [2l). However, in a multiple gateway scenario, 
the precoding scheme shall be reconsidered. 
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Fig. 2. The picture depicts the multiple gateway multicast multiheam satellite structure. In contrast to the single gateway 
architecture, several gateways transmit the data to be delivered to the coverage area. This entails two main drawbacks. First, 
certain inter-gateway link (dotted lines) shall exist. Second, each ground unit must compute an independent submatrix of the 
overall precoding matrix. 


In order to keep the payload eomplexity low, eaeh feeder link reeeiver at the payload will only 
route signals a set of feeds Ng iox g = 1,..., (7. Otherwise, a very eomplex analog seheme shall 
be implemented. In addition, it is important to remark that the feeder link bandwidth resourees 
ean only support a small number of Ng feed signals. Considering this approaeh, the preeoding 
matrix W must be partitioned into G sub-matriees leading to a large deerease of the aehievable 
throughput. In the following subseetions we propose some teehniques in order to overeome this 
limitation. 

In addition, eonsidering that in order to reduee the inter-gateway eommunieation, eaeh gateway 
individually eomputes its preeoding sub-matrix, eertain eooperative seheme shall be eoneeived. 
Indeed, eaeh gateway only have aeeess to the feedbaek from its eorresponding set of feeds 
Ng. With this, eertain CSI eooperation among gateways shall be established so that inter-beam 
interferenee is mitigated. 

A. Precoding Scheme 

As diseussed in the previous seetion, the g-\h gateway only has aeeess to a set of Ng of the 
overall feed elements loeated at the payload. Additionally, the gateway serves a set of Kg beams 
out of the K, leading to a total amount of served users equal to KgQ. We will assume a known 
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feed alloeation per gateway and fixed Ng and Kg. Preeisely, we will eonsider that the feeds are 
assigned in a eonseeutive fashion over the ehannel matrix as it is defined in the following. 

Let us assume that eaeh gateway has aeeess to the overall CSI. With this, the matrix to be 
transmitted through the feeder link of the g-\h gateway jointly with the user symbols is 

G (52) 


In other words, as eaeh gateway only have aeeess to a eertain set of feed elements the symbols 
to be transmitted to the Kg beams shall be linearly transformed with in order to inerease 
the overall throughput. Whenever eaeh gateway has aeeess to the overall ehannel matrix H, 
the preeoding seheme ean be eomputed as deseribed in Seetion III and adapted to the multiple 
gateway seenario. This ean be done by beams of setting zero entries in the preeoding matrix 
whenever the gateway does not have aeeess to a eertain feed. This transformation leads to a 
bloek diagonal preeoding matrix as follows 

0 ■ ■ ■ 0 
0 W2 0 


Wm-gw — 


e C 


NxK 


(53) 


0 0 • • ■ 

Evidently, robust designs ean be applied without any additional penalty. 


B. CSIT Sharing 

Sinee the gateway only has aeeess to eertain beams, its available CSI is redueed Preeisely, 
the g-\h reeeives from its feedbaek link the following matrix 

= H {{{{g - l)QKg + 1) : gQKg, 1 : iV)) G (54) 

where the Matlab notation has been used for the sake of elarity. However, in order to eompute 
the preeoding matrix, the gateways need the ehannel effeet between their assigned feeds to the 
non-intended users. With this, every gateway must transmit over the inter-gateway link (e.g. a 
fiber optie) the information fed baek from its users eorresponding the effeet of the non-assigned 
feeds. Mathematieally, the g-\h gateway must share 


(55) 
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where 

= H ((((^ - l)QKg + 1) : gQKg) , ((/ - l)iV, + 1) : gN^) . (56) 

This eooperation among gateways require a total amount of 

{G-l)QKgNg (57) 

eomplex numbers to be shared by each gateway. As a consequence, it is essential to reduce this 
communication overhead in order to reduce the overall system cost. 

One approach is to limit the sharing between the different gateways and only consider the C 
closer gateways. With this, the overall data overhead reduces to 

CQKgNg. (58) 

Another alternative is to apply certain compression to the transmit channel submatrices. For 
instance, we could use the eigenvector associated to the largest eigenvalue of each matrix H®*. 
This will lead to a total communication overhead of 

(G-l)iV,. (59) 

Evidently, whenever each gateway has a more precise version of the channel matrix, the larger 
throughput can be obtained. In the simulation section this is evaluated and different cooperation 
schemes are evaluated. 


VII. Simulations 

Considering a reference scenario of a geostationary satellite with N = K = 245, we evaluate 
the proposed precoding designs considering a full frequency reuse scenario. Array fed radiation 
pattern has been provided by the European Space Agency and it takes into account the different 
user locations over the European continent. The link budget parameters are described in Table 
1 . 

All results have been obtained considering 500 channel realizations and a phase effect between 
feeds X = 10 degrees. Moreover, throughput values are obtained by means of the user SINR and 
the efficiency (bit/symbol) given a minimum Packet Error Rate (PER) of 10“® considering DVB- 
S2X. It is important to remark that this relationship has been obtained from ll29ll considering 
the PER curves. 
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TABLE I 

USER LINK SIMULATION PARAMETERS 


Parameter 

Value 

Satellite height 

35786 km (geostationary) 

Satellite longitude, latitude 

10° East,0° 

Earth radius 

6378.137 Km 

Feed radiation pattern 

Provided hyESA 

Number of feeds N 

245 

Beamforming matrix B 

Provided by ESA 

Number of beams 

100 

User location distribution 

Uniformly distributed 

Carrier frequency 

20 GHz (Ka band) 

Total bandwidth 

500 MHz 

Roll-off factor 

0.25 

User antenna gain 

41.7 dBi 

G/T 

17.68 dB/K 


The outline of the subsequent subseetions is as follows. First, we show the performanee gain 
of the proposed preeoding sehemes considering perfect CSI and single gateway architecture. 
Second, it is shown that larger throughputs can be obtained if user grouping techniques are 
applied. Third, the impact of imperfect CSI is shown and the convenience of robust designs 
is presented. Finally, the multiple gateway architecture is evaluated so as the proposed inter¬ 
gateway cooperation techniques. Remarkably, for a best practice we also consider a reference 
scenario that consists in 4-colouring scheme where the interference is reduced so as the available 
bandwidth. 

Figure[3]presents the system throughput considering the proposed precoding schemes in section 
III. Both of them, R-ZF and MBIM are compared to the average MMSE design presented in 
lHH, ifTTll . It can be observed that both proposals behave better than average MMSE scheme. 
Specially, MBIM offers larger throughputs than R-ZE over the different transmit power values. 
Moreover, the conventional 4-coloring scheme has the lowest achievable rate due to the bandwidth 
limitations. 
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Average MMSE-Total power constraint 

- ^ - MBIM-Total power constraint 

- ★ - R-ZF-Total power constraint 
I MBIM-Per antenna power constraint 
0 R-ZF-Per antenna power consttaint 
0 Average MMSE-Per antenrra power constraint 

- A - Ref.scenario:4 colors 






A - A -A- 4 - A -4- i-i ^ ' 


[dBW] 



Fig. 3. Per beam average throughput with Q — 2,3 users per beam, perfect CSI, no user grouping. 

A. User Grouping 

As discussed before, whenever the users within the same beam have eollinear ehannel veetors, 
large rates ean be obtained. This is shown in Figure IH Similar to the study in @1, in eaeh beam 
Q = 200 users are uniformly distributed, resulting in an average user density of 0.023 uscrs/km"^ 
inside the 3 dB eoverage edge of every beam. 

It results that in all eases, user grouping inereases the throughput. Speoifieally, for Q = 3, 
the user grouping gain beeomes larger than in the Q = 2 ease. Considering these results and 
the assoeiated low eomplexity operation, it results eonvenient to implement this teehnique in 
satellite systems. 

B. Robust Design 

This subseetion eonsiders the effeet of imperfeet CSI at the transmitter. This is modelled 
with a additive perturbation matrix whose entries are Gaussian distributed with zero mean and 
varianee equal to 1. The perturbation values 7 ^ are eonsidered the same for all submatriees 
{K = 1,... ,K) and in the simulations are obtained eonsidering that eaeh submatrix does no 
exeeed the assumed bound. Additionally, 7 fc is set to 1 and no sensitivity analysis is performed 
due to spaee limitations. Figure [5] shows the performanee of our proposal for different 7 values 
represented in a ratio basis with respeet to the ehannel matrix Frobenious norm. For both oases 
Q = 2,3, the proposed teehnique is able to overoome the imperfeet CSI values at the transmitter 
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Fig. 4. Per beam average throughput with Q — 2,3 users, perfect CSI, user grouping. 

leading large throughputs while keep the preeoding complexity low. 

C. Multiple Gateway Transmission 

In order to evaluate the performance of the proposed multiple gateway schemes, we consider 
the following methods and inter-gateway architectures. We will consider a total number of G = 14 
gateways each of them serving either 17 or 18 beams. 

• Scenario 1 : Individual cluster processing. Each gateway processes its set of beams in¬ 
dependently and only receives the CSI from its corresponding beams. With this, it is not 
possible to mitigate the interference of adjacent beams. This is referred to Individual Cluster 
Processing (ICP). 

. Scenario 2 : Gateway g (respectively for all the gateways) collaborates with 4 gateways 
that serve beams directly adjacent to their beams so that ^f-th gateway receives prefect CSI 
of adjacent beams. This is referred to 4 Gateway Collaboration Processing (4GCP). 

. Senario 3 :Gateway g (respectively for all the gateways) collaborates with 7 gateways that 
serve beams directly adjacent to their beams so that g-\h gateway receives prefect CSI of 
adjacent beams. This is referred to 7 Gateway Collaboration Processing (7GCP). 

. Senario 4 : Gateway g (respectively for all the gateways) collaborates with all gateways 
by means of sharing the singular left vector associated with the largest singular value of 
the gateway channel matrix. This is referred to Maximum SVD of Gateway Collaboration 
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Fig. 5. Per beam average throughput with Q — 2,3 users, imperfect CSI, no user grouping. 


(MSVDGC). 

• Senario 5 : Single gateway seenario so that a unique on ground processing unit is able to 
use all available feeds with the overall channel matrix. This scenario refers to Reference 
scenario (Ref). 

Figure [6] shows the proposed multiple gateway architectures for both the R-ZF and MBIR 
precoding designs. Evidently, as it happened in the single gateway scenario, MBIR provides 
larger overall throughputs than the R-ZF in all cases. In addition, it is observed that the larger 
cooperation is considered, the larger throughputs are obtained. Remarkably, MSVDGC method 
offers a good trade-off between overall throughput and inter-gateway communication overhead. 

VIII. Conclusions 

This paper proposed a two-stage low complex precoding design for multibeam multicast 
satellite systems. While the first stage minimizes the inter-beam interference, the second stage 
enhances the intra-beam SINR. A robust version of the proposed scheme is provided based on 
the novel approach of the first perturbation theory. Additionally, user grouping and multiple 
gateway schemes are presented as essential tools for increasing the throughput in multibeam 
satellite systems. The conceived methods are evaluated in a continental coverage area and they 
result to perform better than the current approaches yet offering a low computational complexity. 
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Fig. 6. Per beam average throughput with Q = 3 users, perfect CSI, no user grouping and different multiple gateway 
architectures with MBIR precoding (Left), R-ZF precoding(Right). 


Appendix A 
Proof of Theorem 1 


After some manipulations, the objeetive function can be written 

trace (w« (^h"hc + !^lj W.. j . (60) 

This objective function can be re-written considering the eigendecomposition of HfcH^ = 
VfcSfcVfc , so that 

l^trace (w"Vr W W.,) . ( 61 ) 

Since the objective function depends on the perturbation matrices and they are unknown by the 
transmitter, the following derivations aim at obtaining an upper-bound of (IMl) . Concretely, an 
upper-bound of both and will be presented. The optimization problem (fTSl) can yield to 
a tractable solution. 

Proceeding with the derivation, the following equality holds 

nf Hfc = Hf Hfc + Af Hfc + Hf Afc + Af A^, (62) 

where it can be observed that the first term is the exact term whereas the rest are perturbation 
terms. In the following, we will consider the perturbation effect of 

AKfc = Af Hfc -f Hf Afc -f Af Afc, 


(63) 
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-H - 


in the preceding design. 

Considering the first order perturbation analysis presented in ETl . the eigenvectors of 
can be written as 

Vfc = Vfc (Rfc + I), (64) 


where 

Rfc = D o (vf AKVf Sfc + SfcVf AK^Vfc ) . 
where the /-entry of D is 


Xf-X,’ 


(65) 


( 66 ) 


for f ^ g and 0 whenever f = g. Xf denotes the /-th eigenvalue of It is important to 

remark that it has been considered that is full rank and, therefore, its null space has 0 

dimension. 

From the problem definition (fTSl) it is possible to bound AK such that 

AKfc ^ + 27fc(Tmax I. (67) 

With this and considering the following lemma: 

Lemma 1. For any semidefinite positive square matrices A, B, C and B ^ C ft holds that 


A o B ^ A o C 


( 68 ) 


Proof: See Theorem 17 of ll^ . ■ 

It is possible to write the following inequality 

Vfc ^ Vfc ^Rfc + , (69) 

where 

R, = efcD o (VfVf Sfc + SfcVfV,) , (70) 

where 

e = fl + 27fc<7max (nf Hfc) (71) 

In case the perturbations of the eigenvectors are consider, the Weyl’s inequality can support 
the approximation 
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Lemma 2 (Weyl’s inequality) Given two Hermitian semidefinite positive matriees Ai and A 2 
with their eorresponding eigenvalues eollapsed in the diagonal matriees Si, S 2 , the following 
inequality holds 

Ssum ^ Si + S2, (72) 

where Ssum is a diagonal matrix whose diagonal entries are the eigenvalues of Ai + A 2 . 

Proof: See ifSTI . ■ 

With this, the eigenvalues of ean be upper-bounded so that 

Sfc ^ Sfc -f efcl. (73) 

With both (l69l) and (TT^ . it is possible to relax the worst-ease maximization in (l39l) . In other 
words, it is possible to use (1^ and (17^ as eigendeeomposition of matrix H^. 

Appendix B 
Proof of Theorem 2 

As for the previous theorem, the main objeetive is to find a bound of the objeetive funetion 
so that the dependenee with respeet to variable disappears. 

Considering the first order perturbation model presented in ifTTI . the required eigendeeompo¬ 
sition of is perturbed by the following matrix 

AQ, = AkWa,. (74) 

This leads to the following approximation of the eigenveetors of 

Zfc = Lfc (Mfc -f I) , (75) 

where are the eigenveetors of Z^ and 

M, = N o (Lf AQ,Lf Tfc + T,Lf AQf L^) (76) 

where is a diagonal matrix eontaining the eigenvalues of Z^ and the g, /-th entry of N is 

-^-, (77) 

f ^kg 

for f f g and 0 whenever f = g. Zk^ denotes the /-th eigenvalue of Z^. 

After some manipulations, we ean lower-bound AQ^ so that 

AQfc + II. (78) 

By means of employing this last inequality and the derivation of the previous appendix, a lower 
bound of the eigenvalues of Zfc ean be established. 



28 


Acknowledgment 

The authors would like to thank Pantelis-Daniel M. Arapoglou for his fruitful comments and 
discussions. 


References 

[1] A. Osseiran, F. Boccardi, V. Braun, K. Kusume, P. Marsch, M. Matemia, O. Queseth, M. Schellmann, H. Schotten, 
H. Taoka, H. Tullberg, M. Uusitalo, B. Timus, and M. Fallgren, “Scenarios for 5G mobile and wireless communications: 
the vision of the METIS project,” Communications Magazine, IEEE, vol. 52, no. 5, pp. 26-35, May 2014. 

[2] B. Evans, “The role of satellites in 5G,” in Advanced Satellite Multimedia Systems Conference and the 13th Signal 
Processing for Space Communications Workshop (ASMS/SPSC), 2014 7th, Sept 2014, pp. 197-202. 

[3] J. Tronc, P. Angeletti, N. Song, M. Haardt, J. Arendt, and G. Gallinaro, “Overview and comparison of on-ground 
and on-board beamforming techniques in mobile satellite service applications,” International Journal of Satellite 
Communications and Networking, vol. 32, no. 4, pp. 291-308, 2014. [Online]. Available: http://dx.doi.org/10.1002/sat.1049 

[4] P. M. Arapoglou, A. Ginesi, S. Cioni, S. Erl, F. Clazzer, S. Andrenacci, and A. Vanelli-Coralli, “DVB-S2x 
Enabled Precoding for High Throughput Satellite Systems,” CoRR, vol. abs/1504.03109, 2015. [Online]. Available: 
http://arxiv.org/abs/1504.03109 

[5] L. Cottatellucci, M. Debbah, E. Casini, R. Rinaldo, R. Mueller, M. Neri, and G. Gallinaro, “Interference mitigation 
techniques for broadband satellite system,” in ICSSC 2006, 24th AlAA International Communications Satellite Systems 
Conference, 11-15 June 2006, San Diego, USA, San Diego, UNITED STATES, 06 2006. [Online]. Available: 
http://www.eurecom.fr/publication/1886 

[6] J. Amau, B. Devillers, C. Mosquera, and A. Perez-Neira, “Performance study of multiuser interference mitigation 
schemes for hybrid broadband multibeam satellite architectures,” EURASIP Journal on Wireless Communications and 
Networking, vol. 2012, no. 1, p. 132, 2012. [Online]. Available: http://jwcn.eurasipjournals.eom/content/2012/l/132 

[7] B. Devillers, A. Perez-Neira, and C. Mosquera, “Joint Linear Precoding and Beamforming for the Eorward Link of Multi- 
Beam Broadband Satellite Systems,” in Global Telecommunications Conference (GLOBECOM 2011), 2011 IEEE, Dec 
2011, PP- 1-6. 

[8] G. Zheng, S. Chatzinotas, and B. Ottersten, “Generic Optimization of Linear Precoding in Multibeam Satellite Systems,” 
Wireless Communications, IEEE Transactions on, vol. 11, no. 6, pp. 2308-2320, June 2012. 

[9] D. Christopoulos, S. Chatzinotas, G. Zheng, J. Grotz, and B. Ottersten, “Linear and nonlinear techniques for multibeam 
joint processing in satellite communications,” EURASIP Journal on Wireless Communications and Networking, vol. 2012, 
no. 1, p. 162, 2012. [Online]. Available: http://jwcn.eurasipjoumals.eom/content/2012/l/162 

[10] V. Joroughi, B. Devillers, M. A. Vazquez, and A. I. Perez-Neira, “Design of an on-board beam generation 
process for the forward link of a multi-beam broadband satellite system,” in 2013 IEEE Global Communications 
Conference, GLOBECOM 2013, Atlanta, GA, USA, December 9-13, 2013, 2013, pp. 2921-2926. [Online]. Available: 
http://ieeexplore.ieee.org/xpl/articleDetails.jsp?amumber=6831518 

[11] A. Gharanjik, B. Shankar, P.-D. Arapoglou, and B. Ottersten, “Multiple Gateway Transmit Diversity in Q/V Band Eeeder 
Links,” Communications, IEEE Transactions on, vol. 63, no. 3, pp. 916-926, March 2015. 


29 


[12] B. Devillers and A. Perez-Neira, “Advanced interference mitigation techniques for the forward link of multi-beam broadband 
satellite systems,” in Signals, Systems and Computers (ASILOMAR), 2011 Conference Record of the Forty Fifth Asilomar 
Conference on, Nov 2011, pp. 1810-1814. 

[13] V. Joroughi, M. A. Vazquez, A. Perez-Neira, and A. Perez-Neira, “Multiple Gateway Precoding with Per Feed Power 
Constraints for Multibeam Satellite Systems,” in European Wireless 2014; 20th European Wireless Conference; Proceedings 
of May 2014, pp. 1-7. 

[14] G. Zheng, S. Chatzinotas, and B. Ottersten, “Multi-gateway cooperation in multibeam satellite systems,” in Personal Indoor 
and Mobile Radio Communications (PlMRCj, 2012 IEEE 23rd International Symposium on. Sept 2012, pp. 1360-1364. 

[15] N. Zorba, M. Realp, and A. 1. Perez-Neira, “An improved partial CSIT random beamforming for multibeam satellite 
systems,” in Signal Processing for Space Communications, 2008. SPSC 2008. 10th International Workshop on, Oct 2008, 

pp. 1-8. 

[16] G. Taricco, “Linear Precoding Methods for Multi-Beam Broadband Satellite Systems,” in European Wireless 2014; 20th 
European Wireless Conference; Proceedings of. May 2014, pp. 1-6. 

[17] D. Christopoulos, S. Chatzinotas, and B. Ottersten, “Frame based precoding in satellite communications: A multicast 
approach,” in Advanced Satellite Multimedia Systems Conference and the 13th Signal Processing for Space Communications 
Workshop (ASMS/SPSC), 2014 7th, Sept 2014, pp. 293-299. 

[18] Y. Silva and A. Klein, “Linear Transmit Beamforming Techniques for the Multigroup Multicast Scenario,” Vehicular 
Technology, IEEE Transactions on, vol. 58, no. 8, pp. 4353^367, Oct 2009. 

[19] D. Christopoulos, S. Chatzinotas, and B. Ottersten, “Multicast Multigroup Precoding and User Scheduling for Frame-Based 
Satellite Communications,” Wireless Communications, IEEE Transactions on, vol. PP, no. 99, pp. 1-1, 2015. 

[20] V. Stankovic and M. Haardt, “Generalized Design of Multi-User MIMO Precoding Matrices,” Wireless Communications, 
IEEE Transactions on, vol. 7, no. 3, pp. 953-961, March 2008. 

[21] J. Liu, X. Liu, and X. Ma, “First-Order Perturbation Analysis of Singular Vectors in Singular Value Decomposition,” Signal 
Processing, IEEE Transactions on, vol. 56, no. 7, pp. 3044—3049, July 2008. 

[22] M. Schneider, C. Hartwanger, and H. Wolf, “Antennas for multiple spot beam satellites,” CEAS Space Journal, vol. 2, 
no. 1-4, pp. 59-66, 2011. [Online]. Available: http://dx.doi.org/10.1007/sl2567-011-0012-z 

[23] M. Aloisio and P. Angeletti, “Multi-Amplifiers Architectures for Power Reconfigurability,” in Vacuum Electronics 
Conference, 2007. IVEC ’07. IEEE International, May 2007, pp. 1-2. 

[24] J. Lizarraga, P. Angeletti, N. Alagha, and M. Aloisio, “Flexibility performance in advanced Ka-band multibeam satellites,” 
in Vacuum Electronics Conference, IEEE International, April 2014, pp. 45^6. 

[25] E. Karipidis, N. Sidiropoulos, and Z.-Q. Luo, “Transmit beamforming to multiple co-channel multicast groups,” in 
Computational Advances in Multi-Sensor Adaptive Processing, 2005 1st IEEE International Workshop on, Dec 2005, 
pp. 109-112. 

[26] A. Silva, Y.C.B.; Klein, “Vehicular technology, ieee transactions on,” in Linear Transmit Beamforming Techniques for the 
Multigroup Multicast Scenario, Oct 2009, pp. 4353^367. 

[27] A. Pascual-Iserte, D. Palomar, A. Perez-Neira, and M. Lagunas, “A robust maximin approach for MIMO communications 
with imperfect channel state information based on convex optimization,” Signal Processing, IEEE Transactions on, vol. 54, 
no. 1, pp. 346-360, Jan 2006. 

[28] D. Palomar, A Unified framework for communications through MIMO channels, Phd Dissertation, 2003. 

[29] “Digital Video Broadcasting (DVB); Second generation framing structure, channel coding and modulation systems for 


30 


Broadcasting, Interactive Services, News Gathering and other broadband satellite applications; Part 2: DVB-S2 Extensions 
(DVB-S2X),” ETSI EN 302 307-2, 2014. 

[30] C. R. Johnson, “Partitioned and Hadamard Product Matrix Inequalities,” JOURNAL OE RESEARCH of the National Bureau 
of Standards, vol. 83, no. 6, pp. 3044-3049, November 1978. 

[31] H. Weyl, “Das asymptotische Verteilungsgesetz der Eigenwerte linearer partieller Differentialgleichungen (mit einer 
Anwendung auf die Theorie der Hohlraumstrahlung),” Mathematische Annalen, vol. 71, no. 4, pp. 441^79, 1912. 
[Online]. Available: http://dx.doi.org/10.1007/BF01456804 


